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The Master Field for the Half-Planar Approximation for Large Matrix 
Models and Boltzmann Field Theory 

I.Ya.Aref'eva and A.P.Zubarev* 

Steklov Mathematical Institute, Vavilov 42, GSP-1, 117966, Moscow, Russia 

In this talk results of study in various dimensions of the Boltzmann master field for a subclass of planar 
diagrams, so called half-planar diagrams, found in the recent work by Accardi, Volovich and one of us (I. A.) are 
presented. 



1. Introduction 

The problem of analytic summation of all pla- 
nar diagrams in realistic models is still unsolved. 
Its solution is closely related with problem of find- 
ing the leading asymptotics in iV x matrix mod- 
els for large A'^ and may have important applica- 
tions to the hadron dynamics In the early 
80-s it was suggested |^] that there exists a mas- 
ter field which dominates in the large A^ limit of 
invariant correlation functions of a matrix field. 

The problem of construction of the master field 
has been discussed in many works, see for ex- 
ample p -pOf . Gopakumar and Gross |Q and 
Douglas |15| ] have constructed the master field 
for an arbitrary matrix model in terms of cor- 
relation functions. There has been a problem of 
construction an operator realization for the mas- 
ter field without knowledge of correlation func- 
tions. Recently this problem has been solved in 
I24I and it was shown that the master fields satisfy 
to standard equations of relativistic field theory 
but fields are quantized according to a new rule. 

In this talk we are going to demonstrate that 
an operator realization for the master field for a 
subset of planar diagrams, so called half-planar 
(HP) diagrams, proposed in gives an ana- 
lytical summations of HP diagrams (see |^ fore 
more details). This construction deals with the 
master field in a modified interaction representa- 
tion in the free (Boltzmannian) Fock space. This 
new interaction representation involves not the 
ordinary exponential function of the interaction 
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but a rational function of the interaction. Corre- 
sponding correlation functions satisfy the Boltz- 
mannian Schwinger-Dyson equations which are 
simpler than the usual Schwinger-Dyson equa- 
tions. In particular in the case of quartic in- 
teraction one has a closed set of equations for 
two- and four-point correlation functions. We 
solve explicitly this system of equations. In the 
case of D-dimensional space-time we get a Bethe- 
Salpiter-like equation for the four-point correla- 
tion function. A special approximation reduces 
this system of integral equations to a linear in- 
tegral equation which was considered p6t in the 
rainbow approximation in the usual field theory. 

A solution of the Boltzmannian Schwinger- 
Dyson equations can be considered as a first non- 
trivial approximation to the planar correlation 
functions. Note in this context that in all previ- 
ous attempts of approximated treatment of pla- 
nar theory were used some non-perturbative ap- 
proximation |jll|,^,^. Topologically diagrams 
representing the perturbative series of Boltzmann 
correlation functions look as half-planar diagrams 
of the usual diagram technique for matrix mod- 
els [Q. We compare numerically the two- and 
four-point HP correlation functions with the cor- 
responding planar correlation functions for the 
one matrix model. For large variety of coupling 
constant the HP approximation reproduces the 
planar approximation with good accuracy. This 
fact gives us an optimism and makes sensitive a 
study of the HP approximation for realistic mod- 
els. 
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2. Half-Planar Approximation for the One 
Matrix Model 

The master field in zero dimensional case is de- 
fined as = a + a+ , where a and a+ satisfy the 
following relation 

aa+ = 1. (2.1) 

This algebra has a realization in the free (or 
Boltzmannian) Fock space generated by the 
vacuum |0), a|0) = 0, and n-particle states \n) = 
(a+nO). 

A free n-point Green's function is defined as 
the vacuum expectation of n-th power of master 
field 



G(°) = (0|</."|0). 



(2.2) 



As it is well-known, the Green's function (2.2) is 



dven by a n-th moment of Wigner's distribution 
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This representation can be also obtained as a so- 
lution of the Schwinger-Dyson equations 

n 

r-(0) _ r^(o) ^(0) 

m— 1 

Interacting Green's functions are defined by the 
formula [El 



Gr, 



(O|0"(l + 5,„i(0))-i|O). 



(2.3) 



In contrast to the ordinary quantum field the- 
ory where one deals with the exponential function 
of an interaction, here we deal with the rational 
function of an interaction. In [ESl it was shown 



that under natural assumptions the form (2.3) 
is unique one which admits Schwinger-Dyson-like 
equations. 

For the case of quartic interaction Sint — 
g(j)^ the Boltzmannian Schwinger-Dyson equa- 
tions have the form 



1=1 



G[Ji^_iGn+k-l-l — g[Gn-kGk+2- 



Gn-k+lGk+1 + Gn-k+2Gk + 



Gn- 



k+dG+k-l] 



The distinguish feature of equations (|]j) is that 
for n > 4 and 2 < k < n—1 the right hand side of 



(2.4) does not contain the Green's functions Gm 
with rn > n. This fact permit us to write down a 
closed set of equations for any G2 and G4 



Go 



1 - 9G2G2 — gGi 



G4 — 2G2 — 2gG2G4. 



(2.5) 



These equations follow from (2.4)for7i = 2,A: = l 
and n = 4, fc = 2. We set Go = 1 to exclude 
vacuum insertions. From ( |2.5[ ) one gets 

(2.6) 



2g^Gl + 3gGl + G2 - 1 = 0. 



Using the Cordano formula for 5 < we have 

G2 = -^Ree*^(Vl - 108.g2 - ieVSg)*- 
v5g 
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1 - 3.g + 16.g2 - lOSg^ + 768/ 



Green's functions for n > 4 are also determined 



from (p. 41). 

Now let us compare the Boltzmann theory 
and the planar approximation for the one-matrix 
model. Green's functions for the one-matrix 
model in the planar approximation are defined 
as 



n2n(g) = lim 



1 1 



DM tr ( M 



,f2n\ 



exp[ 



--tr(M^) - ^tr(M4)l, 



(2.7) 



l=k+l 



where 2^ is a normalization factor. The integra- 
tion in (2/7) is over N x N hermitian matrices. 
According the 't Hooft diagram technique the 
perturbative expansion in the coupling constant 
of the correlation functions (2.7) is represented by 
a sum of all planar double- line graphs . Due to 
the normalization factor 7V~(^+"^ external lines 
corresponding to tr M 2" can be treated as lines 
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of a generalized vertex. We shall call two double- 
line planar graphs topologically equivalent if one 
of them can be transformed into the other by a 
continuous deformation on the plane. 

A planar non-vacuum graph is an HP graph 
if it is topologically equivalent to a graph which 
can be drawn so that all its vertices lie on some 
plane line in the right of the generalized vertex 
tr M 2" and all propagators lie in the upper-half 
plane without overlapping |2^ . Also we shall call 
a planar graph HP-irreducible if it is represented 
as an HP graph in an unique way. A simple anal- 
ysis shows that an HP graph is HP-irreducible if 
it does not contain any tadpole subgraphs. By 
a graph with a tadpole we mean a graph with a 
subgraph which contains two lines coming from 
the same vertex so that after removing of these 
two lines the remaining subgraph becomes dis- 
connected with the rest of the graph. 

In |2^ it has been shown that if one considers 
some special approximation for the planar theory, 
namely so called HP approximation, then Green's 
functions of the one-matrix model in this approx- 
imation coincide with correlation functions in a 
corresponding Boltzmann theory. 
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(2.8) 



In the LHS we select only irreducible HP graphs 
(this approximation we call the HP approxima- 
tion) and in the RHS we omit all graphs with 
tadpoles. Let us sketch a pure combinatorial 
proof of equation (2^) It is based on the 

following two statements. The first one states 
that the planar correlation functions without vac- 
uum insertions in all order of perturbation the- 
ory are represented by a sum of all topologically 
non-equivalent graphs without any combinatorial 
factors. According the second one, all graphs in 
the Boltzmann theory contribute into correlations 
functions (2.3) without any combinatorial factors. 
Now let us compare numerically the HP ap- 



proximation with the planar approximation. The 
explicit formulas for an arbitrary planar Green's 
functions are well known II . On the Table we 
give the results of numerical calculations of the 
HP Green's functions G2, G4 and the planar 
Green's function H4 for the same values of 
the coupling constant g. One can see that the 
answers for HP Green's functions G2, G4 prac- 
tically saturate the planar Green's functions H2, 
H4 in the board interval of the values of g. 

It is instructive to reformulate the Boltzmann 
theory so that it reproduces only tadpole-free HP 
graphs. For this purpose let us defined the fields 
^ and (j) as follows (j) — a + a'^jip — b + b'^, where 
a, a+, b and 6+ satisfy the following relations 

aa+ = 1, bb+ ^ 1, ab+ ^ 0, ba+ = 0. 

This algebra has realization in the Boltzmannian 
Pock space under the vacuum |0) : a|0) = 5|0) = 
0. 

Now let us consider the following Green's func- 
tions 

Fr, = (O|V'0"-2V.(l+5^nt)"'|O), 

Sint = gip : : ip = gip(l)4>tp - gipip. 

The Schwinger-Dyson equations for the Green's 
functions F2 and have the form 

^2 = 1- 9F4 + gF2, F4 = ~gF2F4 + F2. 

From these equations we find 



Fo = 



F4 = 



-1+.9 + v/l + 2.g-3.g2 
2.9(1 - 9) 

1 + g- v/l + 2g-3g2 
29' 



Figure 1 




3. Boltzmann Correlation Functions for D- 
Dimensional Space-Time 

In this section we derive the Schwinger-Dyson 
equations for Boltzmann correlation functions in 
D-dimensional Euclidean space. To avoid prob- 
lems with tadpoles let us consider the two-field 
formulation. We adopt the following notations. 
Let ip{x) = ij+{x) + ij-{x), (l){x) =(j)+ {x)+(t)- {x) 
be the Bolzmann fields with creation and annihi- 
lation operators satisfying the relations 

■0"(a;)V'"*"(2/) = = D{x,y), 

^-{x)4>+{y)^4>-{x)^+{y)^Q. 

whereD{x,y) = J + m^)-^e''''^''~y^ is D- 

dimensional Euclidean propagator. The n-point 
Green's function is defined by 



Fnixi,...,Xn) = {0\lp{xi)4>{x2): 



(3.9) 



Let us write down the Schwinger-Dyson equa- 
tions for the two- and four-point correlation func- 
tions. We have 

(- A +m^)^F2{x, y) = gD{x, x)F2{x, y) 

- gF4{x,x,x,y) + 6{x - y), (3.10) 
(- A +m^)yF4{x, y, z, t) = 

- gFi{y, y, z, t)F2ix, y) + 5{y - z)F2{x, t). (3.11) 
Here we also assume that all vacuum insertions 



are dropped out. We see that equation (3.11 



As a consequence we have a closed set of equa- 
tions which are enough to find F2 and F4. We 
define an one-particle irreducible (IPI) 4-point 
fimction r4(a;, y, z, t) as 

Ti{x,y,z,t)^ J dx'dy'dz'dt'F^^{x,x') (3.12) 

D'^{y,y')D-^{z,z')F^\t,t')T4x',y\z',t'), 
where J-4 is a connected part of F4 

Fi{x,y,z,t) = Ti(x,y,z,t) + F2{x,t)D{y, z). 

Note that in the contrast to the usual case in 



the RHS of (3.12) we multiply JF4 only on two full 
2-point Green functions while in the usual case 
to get an IPI Green function one multiplies an n- 
point Green function on n full 2-point functions. 
From ( |3.10| ) and ( |3.12| ) we have 



r4(p, k,r) = -g-g J dk'F2{p + k- k')x 
D{k')Ti{p + k~k\k' ,r) 



(3.13) 

Equation ( |3.13 ) is the Bethe-Salpeter-like equ- 
ation with the kernel which contains an unknown 
function i^2- As in the usual case we can write 
down F2 in term of the self-energy function S2 

F2^ i 

+ m? + Y,2 

and write equation ( ^.10 ) as an equation for E2, 



does not contain six-point correlation functions. 



'^2{p)=g J dkdqF2{k)D{q)x 

D{p~k~q)Ti{p,k,q). (3.14) 

Equation ( |3.14[ ) is similar to the usual relation 
between the self-energy function S2 and the 4- 
point vertex function for (p^ field theory, mean- 
while equation ( ^.13 ) is specific for the Boltz- 
mann field theory. Equations (3.13) and (3.14) 
are drawn on Fig.l. 
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Figure 2 





a) 




b) 



Equations for r4 and S2 contain divergen- 
ces. To remove them we apply to Boltzmannian 
Green's functions (3^) the standard renormahza- 
tion procedure based on i?-operation. The per- 
turbative expansion of the Boltzmannian Green's 



functions (3.9) is represented by the sum of HP 
Feynman graphs. We draw (/)-propagator by thin 
lines and ■(/'-propagator by thick lines. Only two- 
and four-point graphs may be divergent. Any 
two-point subgraph of the HP graph is also from 
the set of HP graphs and performing contractions 
of the two-point subgraphs of the HP graph one 
gets again an HP graph. But four-point sub- 
graphs of an HP graph may be of HP type and 
may be not. We refer to the later case as to the 
case of H-subgraphs. Examples of divergent HP- 
subgraphs and H-subgraphs are drawn on Fig. 2a 
and 2b, respectively. 

A detailed consideration shows that all di- 
vergent parts can be collected to counterterms. 
Four-point divergent parts of H-type require a 
counterterm : ipijj : il)~ip~ . 

Let us assume dimensional regularization and 
the minimal subtraction scheme with a scale fj,. 
For renormalized correlation functions we have 

Fn{xi, Xn;M{fi), g{p), A(/i), = 
(0 1 V' ) (/)(a;2 ) • • ■ ^(a;„- 1 (a;„ ) X 

where 



fJ-'^Sgijj : (fxf) : ip + fi'^SX : ipip : ip~ip~ 

and e = 4 — D. Here we add a new term to the 
interaction. This A-term modifies only the mass 
term in equation ( S.lOj ) (see for more details). 

The renormalized correlations functions satisfy 
the following renormalization group equation 



d 



d_ 



d 



91nM2 



-7)x 



F„(xi, ...,a;„; M{^),g{p), \iiJ,),fJ-) = 0, (3.15) 
where 



IM 



7 = 



The difference between ( 3.15| ) and the usual re- 
normalization group equation is that the ano- 



malous dimension 7 in (3.15) is not multiplied 
on n/2. There are also a differences in the ex- 
pressions for the beta-functions in terms of the 
counterterms 6g, 6X and Z-^ — 1. In dimensional 
regularization the counterterms are poles in s, 

aj(.g, A) 



Sg 



i=l 



6X 



2-^ pi ■ 



= 1 



E 

i=l 



c»(g, A) 



(3.16) 



As in the usual case one can write the low-order 
expression for the /3-function in terms of ci and 
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ai. But now there is a new expression for go in 
terms of Sg and Z^. We have 

50 ^ fJ.%g + Sg)Z,'^, 



Xo = ^i'{X + SX)Z-^ (3.17) 
Applying /i^ to both sides of ( 3.17 ) one finds 



/3a (A) = (1- 



d 



ding 
d 

9 In A 



)(3Ci(.9,A) - ai(g, A)). 



)(2Aci(.9,A)~5i(.9,A)). 



Taking into account the exphcit form of ci{g, A), 
ai(g,A) and 6i(.g,A) we get 

M9A) = £, + ^£^,^0{g\ X% 



f3,{g,X)^^ + 0{g\ A^). 



(3.18) 



Note that there is only a numerical difference of 
the HP /3 function with the standard f3 function. 
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